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Abstract _ 

A compressible near-wall two-equation model is derived 
by relaxing the assumption of dynamical field similarity between 
compressible and incompressible flows. This requires 
justifications for extending the incompressible models to 
compressible flows and the formulation of the turbulent kinetic 
energy equation in a form similar to its incompressible 
counterpart As a result, the compressible dissipation function 
has to be split into a solenoidal part, which is not sensitive to 
changes of compressibility indicators, and a dilatational pan, 
which is directly affected by these changes. This approach 
isolates terms with explicit dependence on compressibility so 
that they can be modeled accordingly. An equation that governs 
the transport of the solenoidal dissipation rate with additional 
terms that are explicitly dependent on compressibility effects is 
derived similarly. A model with an explicit dependence on the 
turbulent Mach number is proposed for the dilatational 
dissipation rate. Thus formulated, all near-wall incompressible 
flow models could be expressed in terms of the solenoidal 
dissipation rate and scraight-forwardly extended to compressible 
flows. Therefore, the incompressible equations are recovered 
correctly in the limit of constant density. The two-equation 
model and the assumption of constant turbulent Prandtl number 
are used to calculate compressible boundary layers on a flat plate 
with different wall thermal boundary conditions and free-scream 
Mach numbers. The calculated results, including the near-wall 
distributions of turbulence statistics and their limiting behavior, 
are in good agreement with measurements. In particular, the 
near-wall asymptotic properties are found to be consistent with 
incompressible behavior; thus suggesting that turbulent flows in 
the viscous sublayer are not much affected by compressibility 
effects. 

1. Introduction 

Density variation in a turbulent flow can come from 
different sources. Some of these are: (i) isothermal mixing of 
gases of different density, (ii) strong temperature gradient in a 
homogeneous fluid, (iii) reactive flows and (iv) compressibility 
effects in high speed flows. Each of these sources gives rise to 
specific aspects that require modeling if the governing equations 
are to be solved. This study makes an attempt to address the last 
source; that is, the modeling of high speed compressible 
turbulent flows. 

Most studies on compressible turbulent flow modeling 
[1-9] invoke the Morkovin postulate [10] to justify the direct 
extension of the incompressible models to compressible flows. 
The postulate was formulated based on early experiments on 
compressible boundary layers along adiabatic walls and 
compressible wakes, and essentially suggested that the 
dynamical field in a compressible flow behaves like an 
incompressible one. This postulate was used by numerous 
researchers to assure that compressibility effects can be 
accounted for directly by the variable mean density in the 
governing equations alone. In other words, the influences of 
fluctuating density on turbulence mixing are essentially assumed 
to be negligible. The validity and extent of Morkovin's postulate 
were reviewed by Bradshaw [2] and he noted that the postulate 
is appropriate for flows where density fluctuations are moderate. 
Therefore, the postulate is not valid for hypersonic boundary 


Therefore, the postulate is not valid for hypersonic boundary * / 
layers, where the Mach number is five or greater, and for flows/ ^ 
with strong pressure gradient effects, such as shock-turbulent- 
boundary-layer interactions. The latter point was confirmed by 
the studies of Wilcox and Alber [1] and Bradshaw [1 1] and led 
to proposals to have the effects of pressure-dilatation correlation 
modeled in the governing equations [12]. A more recent study 
where density fluctuations are also considered has been given by 
Speziale and Sarkar [13]. Besides these modifications, all 
turbulent compressible flow modeling rely on incompressible 
models. 

Two sources of difficulties arise when incompressible 
turbulence models are extended to compressible flows. One is 
due to compressibility itself and another is associated with the 
turbulence phenomena. In compressible flows, the governing 
equations are coupled and temperature cannot be considered as a 
passive scalar. As a result, all other thermodynamic variables 
adopt new roles. Therefore, mathematically, compressible 
flows cannot be considered as straightforward extension of 
incompressible flows. Furthermore, pressure is only a force 
term in incompressible flows and all disturbances propagate at 
infinite speed. On the other hand, pressure also supports finite 
velocity propagation of disturbances in compressible flows. 

Other complications come from the variable mean density, which 
contributes to increased non-linearity of the governing 
equations, and the fluctuating density, which causes the closure 
problem to become more difficult. 

The second source of difficulties has to do with 
turbulence mixing. Here, even for incompressible flows, many 
problems remain to be resolved [14-17], especially when the 
flow is unsteady and/or three-dimensional [18]. However, 
among the many problems associated with turbulence modeling, 
one stands out as most fundamental and urgently needs 
attention. This is the treatment of the near-wall flow [17]. 
Conventional approach is to invoke the wall function 
assumptions; thus implying that near-wall turbulence is in local 
equilibrium. Even for simple wall shear flows, the assumption 
is riot quite valid because near-wall turbulence is not in local 
equilibrium. Consequently, a low-Reynolds-number treatment 
is necessary in order to obtain results that agTee with 
measurements in the near-wall region [17, 19-21]. The need for 
near-wall treatment of flows with heat and mass transfer has also 
been pointed out [22-25]. This problem is expected to be more 
acute in compressible flow modeling [13] where the non- 
linearity of the governing equations are further compounded by 
the variable mean density. 

The present objective is to model near-wall compressible 
turbulent flows where the coupling between velocity and 
temperature cannot be ignored. As a first attempt, only the 
modeling of the turbulent kinetic energy and its dissipation-rate 
equations is considered. With the assumption of gradient 
transport, the two-equation model could be used to effect closure 
of the mean flow equations. Since the transport equations for 
the heat fluxes and the temperature dissipation rate are not 
modeled and solved, a constant turbulent Prandtl number is 
invoked to relate the heat fluxes to the momentum fluxes. In 
view of this assumption, the present approach only addresses 
the issue of compressibility effects on turbulent mixing and not 
on heat transfer and its interaction with turbulence. An attempt 


Present address: CFD Research Corporation, 3325 Triana Blvd., Huntsville, AL 35805 


on this latter problem will be made after the present model has 
been validated 


2. Proposed Modeling Approach 

With the availability of near- wall models for temperature 
variance and its dissipation rate [24], heat fluxes [25], 
Reynolds-stresses [26] and the dissipation rate of the turbulent 
kinetic energy [27], the time is now ripe for their extension to 
compressible flows. In order to consider the effects of variable 
mean density and its fluctuation on turbulence mixing, it is 
necessary to analyse the exact equations and propose appropriate 
models to effect closure. Two approaches are available. One is 
to propose totally new models for the terms in the compressible 
equations, while another is to attempt to extend the 
incompressible models to compressible flows in a credible way. 
Both approaches involve assumptions that could or could not be 
verified experimentally. Since the present knowledge of 
incompressible flow modeling is quite mature, as a first attempt, 
it is expedient to extend these models to compressible flows. 
This can be accomplished by recasting the compressible 
equations in forms similar to their incompressible counterparts 
so that terms with explicit dependence on compressibility effects 
can be isolated separately, and the incompressible limit can be 
recovered in a straight forward and correct manner. 

Since the turbulent kinetic energy equation or k-equation 
is obtained by contracting the Reynolds-stress equations, this 
means that the recasting of the Reynolds-stress equations should 
be attempted first. In other words, the viscous diffusion and 
dissipation terms in the Reynolds-stress equations have to be 
similarly defined as their incompressible counterpans. This 
suggests splitting the viscous dissipation function into a 
solenoidal pan, which is not sensitive to changes of 
compressibility indicators, and a dilatational pan, which is 
directly affected by these changes 28 . When the Reynolds-stress 
equations are written in this form, three additional terms that 
depend explicitly on compressibility effects are present. The k- 
equation is then obtained by contracting the Reynolds-stress 
equations and its incompressible counterpan is recovered 
correctly when density becomes constant and the additional 
terms vanish identically. An equation that governs the transpon 

of the solenoidal dissipation rate (e) of the turbulent kinetic 
energy (k) is derived and modeled along the line suggested 
above. Again, additional terms that depend explicitly on 
compressibility effects appear in the equation. This equation 
also reduces correctly to its modeled incompressible counterpan 
because the additional terms vanish for constant density flows. 
All models proposed for the k and e equations are expressed in 
terms of this solenoidal dissipation rate. A model with explicit 
dependence on the turbulent Mach number proposed by Sarkar 
et al. [28] for the dilatational dissipation is adopted. Thus 
formulated, the two-equation model is valid for compressible 
flows and approaches its incompressible limit in a straight 
forward and correct manner. 

The systematic approach described above, if proven 
successful, could be used to extend incompressible near-wall 
models for heat-fluxes, temperature variance and its dissipation 
rate to compressible flows. A set of equations governing the 
transport of incompressible heat fluxes has been proposed and 
validated against simple flows with heat transfer [25], while a 
similar set of equations for the temperature variance and its 
dissipation rate [24] has also been validated against boundary- 
layer flows. This means that near-wall heat transfer models 
could also be extended to compressible flows using the approach 
proposed above. However, before this extension is undertaken, 
the asymptotic consistency of these models has to be verified. 
Until such time, the assumption of a constant turbulent Prandtl 
number for near-wall compressible flow is inevitable. 

In the following, the compressible equations are first 


derived, then the near-wall modeling of the k and e equations are 
discussed. In section 6, the two-equation model is used to 
calculate compressible boundary layers on a flat plate assuming a 
constant turbulent Prandtl number. Comparisons with 
measurements [9, 29-31] and other calculations, such as those 
obtained using the k-o> model of Wilcox [8], are carried out to 
assess the importance of density fluctuations on the calculated 
results and, hence, the validity and extent of Morkovin's 
hypothesis. 


d Mean Fl gw Equations 


The compressible mean flow equations are obtained by 
applying Favre averaging to the instantaneous Navier- Stokes 
equations which for Newtonian fluids can be written as: 

dp d f , A 

dr + a ^ (pUi,=0 ’ (1) 
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where tjj = 



2,, 3u k x 
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Uj is the i lh component of the velocity vector, Xj is the i th 
component of the coordinates and p, T, p, p, k, C p are 
pressure, temperature, density, viscosity, thermal conductivity 
and specific heat at constant pressure, respectively. Favre 
decomposition is applied to all variables except p and p where 
conventional Reynolds decomposition is assumed. In other 
words 


Uj = <Uj) + Uj . 

(5) 

t = <e> + e” , 

(6) 

p = p +p' . 

(7) 

p=J + p' . 

(8) 


where u and are the Favre fluctuations and p' and p' are the 
Reynolds fluctuations. If < > is used to denote Favre-averaged 
quantities and the overbar the Reynolds-averaged quantities, 
then the mean equations for compres sible flows c an be obtained 
as follows. The above decompositions (5 )-(8 ) are substitute d 
into (l)-(4) and the resultant equations are averaged over time. 
If the turbulent flow is further assumed to be stationary and the 
mean momentum eauapqn and the Reynolds-stress and turbulent 
kinetic energy, k = “<u- u^, equations to be derived later are used 
to simplify the thermal* energy equation, the turbulent mean flow 
equations become 


J^(p<Uj>) = o 


(9) 
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Symbolically, the above equation can be written as 
Qj = ^ij + ^ ij ‘ P + ^ij + ^ij + ®ij + T« 


( 13 ) 



In these equations, \i ■ \i , K ■ K andCp=Cp have been 
substituted and the mean and fluctuating stresses are given by 


With the exception of Gy and Ty, (13) is similar to its 
njcompressible counterpart P6], For an incompressible flow, 
u. sO, and Gy * Ty = 0. Even under this condition, (13) fails 
to reduce properly to the incompressible equation given in Ref. 
26. _Jhe reason lies in the grouping of the tenns 
(Dy v - p£y + C>y| In order to achieve this incompressible limit 

correctly, a re-arranging of the terms in (d*^ - pe*; + d>y) is 
necessary. If viscous diffusion and dissipation in compressible 
flows are again defined similarly to their incompressible 
counterparts, or 




- 3(Ui> , d<Uj> 
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The quantity, C p (0) + kUfcXUfc) + is the mean total enthalpy 
(H). Thus written, (y)-(ll) reduce to their incompressible 
counterparts exactly when density becomes constant. 


An order-of-magnitude analysis is carried out on (9) - 
(11). The result shows that the underlined terms are of smaller 
order and, as a first approximation, could be neglected compared 
to the terms retained. Thus formulated, the compressible 
equations are identical to the incompressible equations and the 
additional unknowns are the turbulent momentum and heat 
fluxes, just as in the incompressible case. The present approach 
proposes to close these equations assuming gradient transport 

As a first attempt, a near- wall two-equation k-e model is used to 
determine the turbulent viscosity and a constant turbulent Prandtl 
number is invoked to relate turbulent momentum and heat 
fluxes. Therefore, the present model cannot fully account for 
the effects of density fluctuation on turbulent heat transfer. 


4, Modeling of the Turb ulent Kinetic E nergy Equation 


The Favre-averaged transport equation for the Reynolds 
stresses p~(UjUj) could be similarly derived as in the 

incompressible case [16]. That is, the ith fluctuating velocity 
equation is obtained by subtracting the mean momentum 
equation from the instantaneous equation. Repeat the same 
procedure to obtain the jth fluctuating velocity equation. The ith 
fluctuating velocity equation is then multiplied by the jth 
fluctuation velocity and vice versa and the two equations are 
added together and averaged over time. Omitting all the algebra, 
the final exact equation is: 
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then the terms (d** - pey + <X>y) can be re-arranged to give 
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(16) 
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Note that (16) reduces to its incompressible counterpan exactly 
when constant fluid properties are assumed. For compressible 

flows, an extra term pey appears in ( 1 6). In addition, three 

additional terms are founa in <X>*. The term p Ey is a dilatational 
term and could be interpreted as compressible or dilatational 
dissipation. This term is only important for compressible flows. 

It should be pointed out that <I>y is given by (17b) and, as 
a result of this particular partitioning, there are several extra 
terms resulted from compressibility and variable viscosity. 
However, at high Reynolds number, dimensional arguments 
reveal that these extra contributions are notj. important. 1/ 
pressure diffusion is further neglected, then Dy, p £y and d>y 
would assume the same form as their incompressible 
counterparts. Therefore, the high-Reynolds-number 
incompressible models proposed for these terms [16, 32] could 
be straight-forwardly extended to compressible flows. 
However, a model for the compressible dissipation term p £y is 



required to complete closure. For high-Reynolds-numoer 
flows, this compressible dissipation could be assumed to be 
isotropic. As a result, the following model is proposed: 


e u = f- 5 u £C 


where e c = ¥r\ 

3 \dx k ; 


(18) 


The modeling of £ c has been attempted by Sarkar et al. 
[28]. They are the first to realize that the contribution of the 
dilatation al dissipation term is important for supersonic shear 
flows. A simple algebraic model, which is based on an 
asymptotic analysis and a direct numerical simulation of the 
simplified governing equations, has been proposed for E c . Their 
proposal could be modified to become 


e c = ctjM^E , 


(19) 


_ — —i ou. 

where is a model constant, M l = 2k fc 2 , pe * H is the 



dissipation of k and c is the local mean speed of sound. 
Therefore, M t is the local turbulent Mach number. It should be 

pointed out that Sarkar et aL's [28] definition of e° is four times 
larger than the definition given in (19) as a result of a different 
splitting of the terms in (16). Consequently, ot] should take on 
a value equal to 1/4 of that suggested in Ref. 28. Based on an 
analysis of decay of compressible isotropic turbulence, Sarkar et 
al. [28] suggested a value of one for their constant. In other 

words, dj = 0.25. If a\ is evaluated based on compressible 
shear flows, its value would be 0.15. The present study adopts 
<Xj = 0.15 for the analyses of boundary-layer flows. 

The k-equation is obtained by contracting (12) and 
making use of (16) and (18) to simplify the resulting equation 
which can be written as: 


5g^ = D^ + D J+5 p ii + ^-5peii-fP e ii 


l T 

+ 2 G ii + 2 1 ii ’ 


(20) 


It can be seen that the terms, dL C>ij and p e h • 
coefficient, u , appearing in G;; and T « ^ - 0 the 

Furthermore, when p is assumed to be constant and u | -0, the 
last three terms in (20) are identically zero and the 
incompressible equation is recovered exactly. The .modeling of 
TV O" and p* £y could be accomplished by drawing parallels 
^th tilei^ncompressibk counterparts [26, 27]. However, this 
requires knowledge of their behavior in the near-wall region. 

The near-wall behavior of (20) can be analysed by 
assuming Tavlor series expansions about the wall for the 
fluctuating quantities. This analysis is similar to the 
incompressible case [26] except that expansions also have to be 
assumed for p and 6 . The proposed expansions are. 


U M = aj y + a2 y 2 + — * 
v” = bj y + 1>2 y 2 + — » 

w H = cj y + cj y 2 + •••» 
9" = d] y + d 2 y 2 + .... 
p' - y + e 2 y 2 + •• 


(21) 


It should be cautioned that, although the velocity and 
temperature expansions are physically correct, the expansion for 
density is an assumption. As pointed out by Bradshaw [11], the 
fluctuating temperature and density could not go to zero 
simultaneously at the wall. Otherwise, it would lead to a zero 
wall p'. In general, temperature fluctuation is assumed to be 
zero at the wall, while p' is not. Here, the assumption is made 
that p / also goes to zero at the wall, however, its value away 
from the wail is finite. Since p' is taken to be essentially zero 
over the whole field in Morkovin's hypothesis [10], the present 
approach could be viewed as a partial relaxation of that 
assumption. Consequently, die proposed model would not be 
valid for all free-stream Mach number and wall thermal 
boundary conditions. Therefore, one of the present objective is 
to analyse the validity and extent of the proposed two-equation 
model. 

For incompressible flows, b j = 0 is obtained by imposing 
the incompressibility condition and becomes a crucial condition 
in near-wall analysis. This important condition holds the key to 
the present extension of the near-wall incompressible models to 
compressible flows. In order to show that ^ indeed vanishes 

under these conditions, the continuity equation for p is first 
derived, or 

^T + J^P’ u k + P'< U k> + P' u kj = 0 - (22) 

Expansions (21) are then substituted into the above equation. If 
(y k ) = 0 at the wall is used, it can be easily verified that, under 
the assumption of (21), bj = 0 is still a valid condition for 
compressible flows, irrespective of the thermal boundary 
condition. Therefore, the assumed p' expansion facilitates the 
modeling of compressible flows, because all terms in (20) have 
similar forms as their incompressible counterpans except the 

extra Ey term which needs to be analyzed. 

Using definition (18) for £■;, it is easily verified that Ey is 
of order y^T The high-Reynolds-number model (19) also has 
similar behavior near a wall. Therefore, it is proposed to extend 
(19) to near-wall flow without modification, while the near-wall 
balance provided by th«; exact Ey is taken into consideration by 
combining it with the <J>y term. As for p Ey, it could be modeled 
by following the arguments presented in Refs. 26 and 27 tor 
incompressible flows. In essence, Refs. 26 and 27 argue that 
the incompressible £y can be set equal to 2 e and the near-wall 
corrections proposed for Ey have little or no effects on the 
behavior of Ey in the region near a wall. This means that Ey can 
be approximated by 2 z in the whole field. In view of this, the 
model for Fey can be assumed to be given by 2p £. Based on 
this model, equation (15) and expansions (21), it can be easily 
shown that the leading term of £ in the near- wall region is a 
constant equal to its wall value E w . Again, the behavior is 
similar to its incompressible counterpart 

Near-wall analysis again shows that turbulent diffusion is a 
higher order term and its high-Reynolds-number model could be 
adopted because it does not affect near- wall balance of the k- 
equation. Consistent with the assumption of gradient transport 
for two-equation models, the incompressible model for turbulenj; 
diffusion of k is extended to compressible flows by writing D k 
= 3((pfi/o k )3k/ax i )/ax i , where Gjjs a constant and \l t is the 
turbulent viscosity defined by jit - ^ ( e * . s 

definition, Cu is a model constant while fy is a damping function 
to be definea later. Based on (21), the leading order term of k in 
the near-wail region is y^. Since £ = ey in this region, k /£ has 
to be of order y*. If the shear stress is defined with respect to 
pT then it can be shown that the lead in g order term of the shear 
stress has to be of order y^ in the near- wall region. Therefore, it 
follows that v\ * ]T/p" is also of order y^ near a wall and this. 



in turn, leads to a similar behavior for the modeled Dj term in 
the near-wall region. This behavior is consistent with the 
behavior of the exact term appearing in (20). In other 
words, the modeled Dr does not affect the near-wall balance of 
(20). k 

According to (18), Ejj = 2e c . As such, the near-wall 
behavior of the exact z[ x is not properly accounted for by the 
proposed model. In the above discussion, it is argued that the 
near-w^ll behavior of could be modeled together with the 
term <!>„. In order to analyse the near-wall behavior of the 

combined teim (<&*\ + p - ^), the behavior of Gjj and Tjj ne* ’ a 
wall has to be studied. The appearance of mean pressure in Gy 
makes the analysis slightly more difficult However, the 
difficulty could be circumvented by making use of the mean 
momentum equation (10). The final analysis shows that the 
combined (Gjj + Tjj) term has the following near-wall behavior, 
namely, 


Gil + T n 0(y 2 ) ; G 33 + T 33 -> 0(y 2 ) ; 

G 22 + T 22 0(y 3 ) . (23) 

Thjs means that to the lowest order, the near-wall behavior of 
(Ojj + pe-) is similar to its incompressible counterpart [26]. 
For incompressible flows, the term, d>jj, can be written into a 
pressure diffusion part and a pressure redistribution part. 
Pressure redistribution is identically zero and since pressure 
diffusion is relatively small, it is usually neglected. Such is not 
the case for compressible flows. The term, C>jj, can again be 
partition into a pressure diffusion pan, which could be 
neglected, and a term involving pressure-velocity-gradient 
correlation. This latter term does not vanish because fluid 
volume changes as a result of density variation. ^Therefore, an 
argument could be made to model the term, (Ojj + p e[j), to 
account for dilatational effects only. In view of this, the 
following model is proposed, or 

(Oj* + pe*) = -ypkj^-j , (24) 


where y is a model constant. 



with P = - 



<e> 

p 


(27) 


where P equals to unity for an ideal gas. 

The near-wall behavior of the modeled k-equation can 
now be analysed using expansions (21). It can be easily shown 
that in the region very near a wall, the modeled k-equation is in 
balance up to order y. Consequently, it does not need further 
modifications to achieve a consistent asymptotic behavior near a 
wall. 


5. Modelling of the Dissip ation- Rate Equation 

The exact transport equation for the solenoidal dissipation 
rate (p'e) can be derived as in the Reynolds- stress equation (12). 
It has been pointed out that the e-equation is the most difficult to 
model even for incompressible flows [13, 16, 17, 26, 27, 32]. 
The reason being that many of the terms in the exact equation arc 
either not known or could not be measured accurately at present. 
Consequendy, the incompressible e-equation is modeled in an ad 
hoc manner to resemble the k-equation in form so that the right 
hand side of the e-equation again consists of four terms; namely, 
viscous diffusion, turbulent diffusion, production and 
destruction of e. The equation is further modified for near- wall 
flows by adding an extra destruction term ^ so that the modeled 
equation remains balance as a wall is approached. There is a 
lack of measurements in compressible flows, therefore, a 
rigorous modeling of the compressible e-equation is not possible 
at present. An alternative is to extend the high-Reynolds- 
number incompressible models to compressible flows and then 
seek a near-wall correction to the modeled e-equation along the 
line suggested in Ref. 27. In view of this, the exact transport 
equation for the dissipation rate is not in a convenient form to 
work with. The proposal of Speziale and Sarkar [13] with the 
dilatational effects explicitly written out will be more 
appropriate. 

Following Speziale and Sarkar [13], the modeled 
transport equation for t with near-wall correction is written in 
the simplified form; namely 


The proposed models still fail to close the k-equation 
because of the presence of in Tjj and Gjj. Therefore, it is 
necessary to shed some light on the modeling of u., which is 
identically zero for incompressible flows. Using Favre 
averaging, it can be shown that -p'u. = p u. In other words, 

m m m 

u. = - p'u. /p. Previous proposals for -p'u. are based on the 
gradient transport assumption; namely, 



(25) 


p— ) + D^ + P E - , (28) 

3xJ J d*i 

where D* is the turbulent transport of £, P e is the production of e 
due to aeviatoric strains, A e is the destruction of z and £, is a 
near-wall correction for compressible flows. The second last 
term on the right hand side of (28) is exact and results from the 
writing of the exact e-equation into the form of (28). When the 
dissipation-rate equation is formulated in this form, it is 
reducible exactly to its incompressible counterpan and, 
therefore, the terms D« P c and Ag can be modeled by a variable 
density extension of their incompressible models. Following the 
suggestions of Refs. 13 and 27, the models proposed for D e , P E 
and Ag are; 


Dpe _ 


where Op is a model constant. However, a more elaborate way 
to model the term is to adopt the proposal, 


u'=- 0 -= Cp=Mu’u*>|jr 

p pe 1 3 ® x j 


(26) 

also 


Dfc- 


d(n,de) 


dxi\ G £dxiJ 


(29a) 


(29b) 


A e = Cc2pS- , 


where Cp is a model constant Alternatively, the term can 
be modeled by 


(29c) 


where the model constants Q] and Q 2 are the^same as those 
given in Ref. 27 for incompressible flows and e = e - E w . It 
should be noted that the mean dilatadonal effects are accounted 
for exactly by (29b) for compressible flows and that these 
models reduce exactly to their incompressible counterparts when 
the flow Mach number becomes very small. In addition, the 
ordering of these model terms is similar to their incompressible 
counterparts. Therefore, the near-wall function \ can be 
determined in a manner similar to that proposed in Ref. 27. 


« ( 1 + 3.45/flQ tanh (yVl 1 5) . (32) 

where y+ = yu x fv is the wall coordinate and u^x) is the friction 
velocity. In this definition, y is taken to be the normal 
coordinate and x the stream coordinate. 

6. Model Validadon 


The incompressible form of (28) with model terms given 
in (29) is identical to that proposed in Refs. 26 and 27. In these 
studies, the coincidence condition suggested by Shima [21] was 
used to determine This is equivalent to requiring the modeled 
e-equation to achieve balance behavior in the near-wall region at 
least up to order y. The approach used to deduce q is to 
assumed a functional form for q with two undetermined model 
constants. One of the constant can be determined from near- wall 
analysis, while the other is evaluated using computer 
optimization. The £ function thus determined has been used in 
Ref. 27 to calculate flat plate boundary-layer flows and in Ref. 

26 to calculate fully-developed channel and pipe flows. These 
calculations were carried out over a wide range of flow 
Reynolds number. The results were compared with direct 
simulation data as well as measurements. Very good agreement 
has been found for both the limiting behavior of the turbulence 
quantities and e when compared to direct simulation data [33- 
35]. Furthermore, the two-equation model calculations of Ref. 

27 arc found to give better results than those obtained by 
Speziale et al. [36]. In view of this success, the same approach 
can be used to determine % for compressible flows. 

The functional form assumed in Ref. 27 is adopted here, 
or 


5 = fw,2P 




(30) 


where f w> 2 is a damping function that goes to one at the wall and 
zero far away from the wall. It is defined in Ref. 27 as f w 2 = 
e -(Rt/64)^ w here R t = k 2 /Ve is the turbulent Reynolds number. 
The function e* is defined as t = e - 2V1C/Y 2 by gencralizin^the 
incompressible definition used in Ref. 27. Similarly, e is 
defined with e w given by e w = 2v (dVF/9xj£\ Once ^ is 
postulated, the near- wall behavior of (28) and the modeled terms 
of (29) can be analysed using expansions (21). If the modeled 
equation is again required to be in balance up to order y, then it 
can be easily shown that N = 2 - C E 2 because the mean 
dilatational terms are of order y. Therefore, to order y° they do 
not contribute to In Ref. 27, the pan involving C e 2 in N is 
grouped together with M to give Mj = + M) and its 

value is determined through computer optimization studies. 
Again, the same procedure is followed in the present study to 
determine Mj. 


A first step to validate the two-equation model for 
compressible flows is to apply it to calculate fiat plate boundary- 
layer flows with different wall boundary conditions and free- 
stream Mach numbers. In this initial attempt, heat flux is not 
modeled separately. Instead, it is related to momentum flux via 
the assumption of a constant turbulent Prandtl number. The 
rationale for doing this is to carefully assess the assumption 
made in (21) concerning the expansion for p\ whose validity 
affects the near-wall analysis used to justify the direct extension 
of the incompressible near-wall function % to compressible 
flows. By choosing the simplest type of compressible flows to 
validate the model, a careful analysis of the validity and extent of 
the expansion for p' can be carried out A second objective of 
this validation is to determine, if possible, the validity and extent 
of Morkovin’s hypothesis. In other words, it is hoped to 
evaluate the Mach number range and the type of wall thermal 
boundary conditions in which the effects of fluctuating density 
can be neglected in the modeling of the governing equations. 
Therefore, the experimental data chosen for comparisons are 
selected from three different groups; one with adiabatic wall 
boundary condition and varying free- stream Mach number, 
another with constant wall temperature and varying frce-stream 
Mach number, and finally the variations of skin friction with 
free-stream Mach number and wall temperature as prescribed by 
the van Driest II formulae given in Ref. 9. The data for the first 
two groups are selected from Ref. 31. Only one case with 
constant wall temperature is selected. The reasons being that the 
other cases reported in Ref. 3 1 are cither not accurate as far as 
the measured skin friction is concerned or the measured mean 
velocity and temperature are doubtful. In the following, the 
governing equations for compressible flat plate boundary- layer 
flows are first presented. Then the calculations and comparisons 
with data are examined and the validity of the two-equation 
model is studied in detail. Finally, the model performance and 
its proposed improvement are discussed. 

6.1 Gove rning Equations 

Two-dimensional, steady compressible boundary layers 
on a flat plate are considered. If the usual boundary-layer 
approximations are made, then equations (9) - (11), (20) and 
(28) can be substantially simplified. For the sake of 
completeness, the boundary-layer equations in Cartesian x-y 
coordinates are listed here as: 

^(p<U>) + ^(p<V>) = 0 ’ (33) 


Finally, to complete closure of the governing equations, 
gradient transport is assumed for the Reynolds stresses and the 
relation is given by 
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In addition, a constant turbulent Prandtl number is assumed so 
that heat and momentum diffusivity can be related by Pr t = v ^fa 
) t where a t is the turbulent thermal diffusivity. The damping 
function f„ appearing in the definition of p t can now be defined. 
In view of the similarity of the present k and e equations with 
those for incompressible flows, the f^ used in Ref. 27 is 
adopted here. It is defined as 
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In writing down these equations, (33) has been used to relate the 
shear stress to the mean velocity gradient It is also noted that 
Gjj is zero for flat plate boundary-laycr_flows while the only 
term of importance in is that given by u" . If (25) is used to 
evaluate u" t then Tjj is of order y 2 and is not important in the 
near-wall region, which according to So et al. [27] is bounded 
by 0 £ y+ £ 5. The model constants are taken from Refs. 27 
and 32 and are given by: Cm = 0.096, Cm = 1.5, C^ ® 1.83, 
a k = 0.75, G e = 1.45, a! « 0.15, a p « 0.5 and y = 0.182. The 
turbulent Prandtl number Pr t is specified as 0.9, while Pr is 
assumed constant and taken to be 0.74. Sutherland and power 
laws are used to evaluate fluid viscosity at an appropriate 


reference temperature. The constants c p and <Xj are determined 
by calibrating the calculated results against some classic 
experiments such as those provided in Refs. 29 and 30. Once 
determined, they are kept constant for all other calculations. In 
fact, the cases to be compared in the next section are different 
from the cases used to calibrate Gq and <Xj. As for N and Mj, 
the final choice of values adopted are 2 and 1.5, respectively, 
just as in Ref. 27. 


The boundary conditions are no slip at the wall for mean 
velocities and k, and zero heat flux or constant wall temperature 
at the wall for enthalpy. As for e, its value at the wall is given 
by 2v (OVlc/Sy) 2 . At the edge of the boundary layer, free-stream 
conditions are w specified for both mean stream velocity and 
enthalpy. In principle, k and e should be zero in the free stream. 
However, in practice, they are assumed to take on some very 
small values, of the order of 1 0' 7 , in the free stream. Thus 
formulated, the above equations and the appropriate boundary 
conditions can be solved numerically using the boundary-layer 
code developed by Anderson and Lewis [30] and modified by 
So et al. [27]. 


All measurements used to validate the near-wall two- 
equation model are drawn from Ref. 3 1 . The calculations are 
carried out over the range, 0 < M^, < 10, for adiabatic wall 
boundary condition and over the range, 0 < 0 w /0 r < 1, for 
constant wall temperature condition. Here, 0 W is the wall 
temperature and 0 r is the recovery temperature for adiabatic wall 
boundary condition and is the adiabatic wall temperature for 
constant wall temperature. It should be pointed out that 0w/0r = 
1.0 corresponds to adiabatic wall boundary condition while 
0 w /0 r < 1 indicates that the wall is cooled. Since only mean 
flow properties are available from Ref. 31, comparisons are 
made with these measurements and another set of calculations 
using the k-co model of Wilcox [8]. All comparisons are made at 
the same momentum thickness Reynolds number (Rq) as the 
measurements. Four sets of data are chosen and these are cases 
55010504, 53011302 and 73050504 with adiabatic wall 
boundary condition and case 59020105 with constant wall 


temperature. The free- stream Mach number (M^) for these 
cases are 2.244, 4.544, 10.31 and 5.29, respectively, while the 
corresponding R e ’s are 20,797, 5,532, 15,074 and 3,939. 
Therefore, the data cover a wide range of Ra and The 
variations of skin friction with and 0 w /0 r are compared 
with the van Driest II formulae reported for a fixed R^ as 
specified in Ref. 9. Finally, an assessment of compressibility 
effects on near-wall flows is attempted by comparing the 
calculated turbulence statistics for the different cases 
investigated. 


6.2 Comparisons with Data 

The results are organized in the following manner for 
presentation. Mean stream velocities are normalized by tu to 
give u+ and thc^are plotted versus In y* (Fig. 1), where y* is 
defined as i^y/v w . A direct plot of the mean velocities is also 
given in Fig. 2 where <U>/U« versus y/5 is shown. Here, 5 is 
the boundary layer thickness defined as the location of y in 
which <U>/U« = 1.0 as specified by the measurements and U„ 
is the free-stream velocity. On the other hand, mean 
temperatures are normalized by 0 W the free-stream temperature, 
and are shown versus y/S (Fig. 3L Plots for thc_properties, k+ 
= k/u^, e + = e v/iLp -uv + = -uv/ur and -0v + = -Gv/U**©*,, are 
presented in-terms of y* and y/5. Near-wall behavior of k+, e + , 
-uv + and -0v + (Figs. 4 - 7) are discussed first, then the 
distributions of Ir 1 * and -uv + in the outer region are examined 
(Figs. 8 and 9). In Figs 8 and 9, 5 is not interpreted from 
measurements; rather it is evaluated at the y location where 
<U>/U oe = 0.9974. Only the budget of k in the near- wall region 
for case 73050504 (M^ = 10.31 and © w /0 r - 1.0) is presented 
(Fig. 10) because the k budgets for the other cases are 
essentially similar to that shown in Fig. 10. The effects of Mach 
number on the asymptotic behavior of k are examined by 
plotting versus Mo* (Fig. 11), where a k is the leading 
coefficient in the expansion of k+ in terms of y^,. According to 
Ref. 27, = 2a*. Therefore, by examining a k versus M«,, the 

variation of £* with Mach number is also evident. Other 
asymptotic properties are tabulated in Table 1 for comparison. 
The variations of skin friction coefficient, Cf = 2i w /p U*,, 
with Moo and ©^/©r are compared with van Driest II results [9] 
in Fig. 12. Finally, the mean velocity plots in terms of the 
compressible u + defined as u* = (j(p/p^)d<U>)/u T for two 
different M« are shown in Fig. 13.® This figure is provided to 
illustrate the deviation or lack thereof from the van Driest law of 
the wall for compressible flows [37, 38]. 


Two versions of the present k-e model are used to 
calculate boundary-layer flows. One designated k-e model/ 1 
solves the k and e transport equations as given in (38) and (39). 
The second designated k-e model/2 solves (38) and (39) with all 
additional compressible terms neglected and the dkjdy term 
omitted in (37). In other words, the two-equation model for k-e 
model/2 is a direct variable density extension of the two-equation 
incompressible model of Ref. 27. These calculations can be 
used to evaluate the validity and extent of Morkovin’s 
hypothesis and the importance of having an asymptotically 
consistent near-wall correction for two-equation models. 

Four sets of u + results are shown in Fig. 1. In the figure, 
the calculated and measured Cf and the Cf determined from the 
van Driest II formula of Ref. 9 are listed for comparison. The 
log-law shown is used to demonstrate the existence of a log 
region in the calculated and measured flows, while the von 
Karman constant K is taken to be 0.41. It is recognized that the 
intercept is a function of Mach number; however, in this figure, 
the intercept is taken to be 4.7. The actual value used is not 
important because the purpose here is to illustrate the slope of 
the log-law. It can be seen that a log region indeed exists for all 
calculated and measured flows. The slope is fairly constant for 
the three adiabatic wall cases tested and the k thus determined is 
approximately 0.41. For the cooled wall case, the k-e model 
predicted slopes are slightly different from that calculated by k-co 
model. None of these slopes yields a von Karman constant of 
0.41 though. Cf is predicted correctly by all three models with a 


maximum error of less than 5%. For the cooled wall case, the 
measured Cf is substantially higher than the van Driest II value 
and, according to Ref. 31, is not as accurate as the measured Cf 
for the other cases studied. The model calculations are in good 
agreement with the van Driest II values for ail cases examined. 

Calculated u + profiles correlate well with measurements. 
At high Mach numbers, there seems to be substantial difference 
between model predictions and measurements in the outer region 
of the boundary layer. Such a difference is also noted when the 
plots are given in terms of <U>/Uoo versus y/5 in Fig. 2. From 
these plots, it can be seen that k-E model predictions of <U>/U« 
are in better agreement with data than k-oo calculations (Fig. 2), 
The agreement between k-E model predictions and measurements 
is good up to = 10.31. On the other hand, the discrepancy 
between calculations and measurements deteriorates as M« 
increases for the k-co model. Therefore, the proposed k-E 
models represent improvements over existing models whose 
predictions are correct only for £ 5 (see e.g. Ref. 39), 
Since most existing models do not have an asymptotically 
consistent near-wall correction and are direct extensions of their 
incompressible counterparts, the discrepancies display by these 

models for > 5 are understandable. Present results show 
that, if the near-wall flow is modeled in an asymptotically cornea 
and consistent manner, the incompressible models can be 
straight-forwardly extended to compressible flows with a free- 
stream Mach number as high as 10, It should be pointed out that 
the k-co model is not an asymptotically consistent near-wall 
model. Therefore, its prediction of <U>AJ W in paticular that 
for the cooled wall case, is not as good as k-E model 
calculations. One reason could be the fairly low R* (3,939) for 
this case. Since the results of k-e model/2 arc also in good 
agreement with measurements (Figs. 1 and 2), the comparisons 
suggest that the additional compressible terms in (38) and (39) 
are not too important. However, an examination of the mean 
temperature results tends not to support this conclusion (Fig. 3). 

The mean temperature profile comparisons are shown in 
Fig. 3. Predictions by the k-co model show substantial 
discrepancies compared to the k-E model calculations and 
measurements; particularly for the cooled wall case (Fig. 2a). 
Discrepancies between measurements and k-co predictions 
increase as M«, increases and as 0w/B r decreases. On the other 
hand, the agreement between k-E model/1 predictions and data 
improves as M* increases for adiabatic wall boundary condition. 
This is not true for k-e model/2 where the disagreement with 
data is quite substantial at M« = 10.31 (Fig. 2b). The 
predictions of the cooled wall case (Fig. 2a) tell a different story. 
It seems that k-e model/2 gives as good a prediction of the 
cooled wall case as that of k-e model/1 (Figs, la and 2a). The 
following three reasons could be put forward to explain this 
behavior. Firstly, the p’ expansion may not be totally valid for 
constant wall temperature boundary condition. Secondly, the 
proposed compressible models may be more applicable for 
adiabatic wall boundary condition. Thirdly, the assumption of a 
constant turbulent Prandtl number may not be appropriate. 

The near-wall distributions of k+ for the four cases are 
shown in Fig. 4. Only the predictions of k-e model/1 and k-co 
model are compared. The calculations of k-E model/2 are^not 
shown; instead, the limiting behavior of k*, -uv+ and -0v + is 
tabulated in Table 1 for comparison. In general, the predictions 
of k-co model are substantially lowered than those of k-E 
model/1. The peaks are about 40% lower than those predicted 
by k-e model/ 1 and the locations where the peaks occur are 
calculated to be further away from the wall than k-e model/ 1 
predictions. According to k-e model/1, the peak of k* decreases 
as Moo increases. The decrease is more than 20% over a Mach 
number range of 10. On the other hand, a slight cooling of the 
wall at Moo = 5.29 causes the peak of k+ to decrease to the same 
level as that for the case of adiabatic wall with Moo = 10.31. 
These results suggest that wall cooling has more influence on 
reducing turbulent mixing compared to compressibility effects. 
Overall, compressibility reduces turbulent mixing and the 


reduction increases with Moo- The near- wall distributions of E + 
as calculated by k-e model/1 are plotted in Fig. 5. It can be seen 
that the distributions are very similar to those shown in Ref. 27 
for incompressible flows. The variations of with Moo and 
© w /© r are very similar to those of k + . Again, maximum e* 
occurs at the wall and a plateau in e + is found in the range, 7 £ 
y+ £ 13. This means that compressibility has little or no effects 
on the near-wall behavior of e + . The exception is that increases 
in compressibility and wall cooling tend to decrease e£. 

In general, k-co model gives a very accurate prediction of 
-v\ y near a wall. Its predictions are as good as those given by 
k-e modcl/1 (Fig. 6). From this set of predictions t jhe following 
observations can be made. Firstly, the peak of -uv+ decreases 
with increasing M« and decreasing 0 w /© r . Secondly, as M^ 
increases and © w /0 r decreases, the location of the peak moves 
towards the wail. Thirdly, the rate of decrease of -uv+ in the 
range, 30 £ y* £ 100, increases as _M« increases. Finally, the 
asymptotic near-wall behavior of -uv+ is listed in Table L-for 
comparison. Much the same behavior is also true for -0v + 
whose distributions in the near-wall region are shown in Fig. 7. 

If the distributions of 0v* are plotted instead of -9v+, the curves 
will have the same shape as those shown for -uvl — Therefore, 
the observations drawn for -uv+ are also valid for 0v+. 

The distributions of k + and -uv + across the boundary 
layer are compared in Figs. 8 and 9, respectively. In all cases 
shown, k-co model over-predicts k+ and -uv+ in the outer part of 
the boundary layer compared to the calculations of k-e model/1. 
The over-prediction extends across the range, 0.2 £ y/5 £ 1.0. 
Reduction of turbulence activities in the outer pan of the 
boundary layer is clearly evident when either compressibility or 
wall cooling effects are present. The reduction increases as 
increases and ©w/0 r decreases. Therefore, it is expected that 
turbulence activities will be substantially reduced in a flow 
where the free-stream Mach number is large and the wall is 
highly cooled. 

The near-wall k budget for case 73050504 is plotted in 
Fig. 10. Other budget plots are not shown because they are 
essentially similar to that given in Fig. 10. It can be seen that the 
k budget bears a lot of similarity with that calculated for 
incompressible flows (see e.g. Ref. 27). The additional 
compressible terms have negligible effect on the near-wall k 
budget. Therefore, the assumptions made to derive the near- 
wall function £ in the dissipation-rate equation are justified. 
Again, viscous diffusion balances dissipation at the wall. This 
balance extends to about y* = 4 where turbulent diffusion and 
production become important. In the region, 4 < y* £ 15, 
viscous and turbulent diffusion, production and dissipation are 
equally important Beyond y* = 15, production and dissipation 
are in balance, just as in the case of incompressible flows. 
Consequently, the near-wall k behavior is very similar for both 
incompressible and compressible flows. 

According to Refs. 25 and 27, Taylor series expansions 
about y+ = 0 can be assumed for k + , -uv+ and -0v+. For 
incompressible flows, the expansions are valid up to about y w * 
7. This range may not be applicable for compressible flows. 
Nevertheless, such expansions for small y* can still be 
assumed. With the help of (21) these expansions can be written 


as: 

k+ = a k (y+) 2 + b lc ( y ;)3 + , (40a) 

= a uv (y*) 3 + b uv (y* ) 4 + (40b) 

-07+ = a ve (y^ + b ve (y+ ) 4 + ...... , (40c) 


[III I I Mill ll Hill I lllillllillln IIInlLi Hill H Iii,| 


where the a's and b's are time-average coefficients that are 
functions of x. A similar expansion can be deduced for e+. 
Again, using (21), the definition of e and its wall boundary 
condition, the expansion for can be written as: 


e + = 2a k + 4b k y; + (41) 


From these expansions, it can be easily deduced that 
k +2 /e + (y*) 2 = 0.5. Therefore, the asymptotic behavior of 
k +2 /e + (y^) 2 is 0.5 and is independent of M« and wall thermal 
boundary conditions. The accuracy in which a model can 
predict this quantity is a reflection of the asymptotic consistency 
of the model. Table 1 shows that k-e model/1 is indeed 
asymptotically consistent while k-e model/2 is not as good. As 

for the k-co model, its prediction of this limiting value is poor, 
therefore, it is not listed in Table 1. 

The "a" coefficients can be determined from the 
calculations and their values are also listed in Table 1 for 
comparison. It can be seen that a^ varies with free -stream Mach 
number. A plot of a^ versus for adiabatic wall boundary 
condition is shown in Fig. 11. The value of a^ for the 
incompressible case is taken from Ref. 27 and is plotted at M« = 
0. Clearly, the trend is to approach an asymptotic value for a* at 
high M«. This decrease in a k is one of the reason why k 
decreases for high Mach number flows (Figs 4 and 8). The 
physical reason is that compressibility tends to hinder turbulence 
mixing. As a result, both turbulent shear stress and kinetic 
energy decrease significantly as M«, increases (Figs. 4, 6, 8 and 
9). Since = 2a* according to (41), dissipation at the wall is 
also dependent on M«. There is no clear trend for a uv and a ve . 
However, the values of are consistent with those calculated 
for incompressible flows [27] and direct simulation data [35]. 
The value of a v0 is essentially zero. Since there is no data 
available, its correctness cannot be verified. 

Finally, the ability of the k-e models to predict skin 
friction coefficient over a range of M« and wall temperatures is 
illustrated in Fig. 12. In Fig. 12a, the variation of Cf/(Cf)j with 
M« for the case of adiabatic wall boundary condition is shown. 
Here, (Cf)j is the skin friction coefficient for an incompressible 
flow evaluated at = 10 4 and is determined to be 2.73 x 10’ 3 . 
The figure shows a comparison of the calculations of k-E 
model/1 and k-e model/2 with the van Driest FT distribution. 
Below Moo = 5, the calculated variations of Cf/(Cf)j with M^ are 
slightly lower than the van Driest II distribution but they are 
slightly higher beyond M« = 5. Essentially, there is no 
difference between the predictions of k-e model/1 and k-e 
model/2. This means that both versions of the k-e model give a 
correct prediction of the Cf/(Cf)j variation with M« for adiabatic 
wall boundary condition. The predictions for the cooled wall 
case are not as good, especially at low temperature ratio (Fig. 
12b). Three sets of calculations are presented. These are k-e 
model/1, k-e model/2 and a third version of k-e model/1 with the 
dk/dy term in (37) neglected. Calculations for this case are 
carried out at M« = 5, R^= 10 4 and the incompressible Cf is 
again determined to be 2.73 x 10 -3 . It can be seen that error of 
5% or larger starts to accumulate at approximately 0 w /0 r = 0.4 
for k-e model/1. This trend is contrary to previous calculations 
[8]. An examination of the governing equations solved by other 
researchers revealed that, besides differences noted in the 
turbulence model equations, the mean energy equation solved by 
these researchers does not include the term dk/dy in the right 
hand side of (37). Indeed, when the dk/dy term is neglected, an 
overall improvement is obtained. The predicted Cf at 0w/9 r = 
0.2 is increased by about 6%, thus giving a better agreement 
with the van Driest II formula. If the additional compressible 
terms in the k-e equations are further neglected (k-e model/2), 
the calculated Cf is only improved by about 3%. The remaining 
disagreement could be attributed to the assumption of a constant 


turbulent Prandtl number. When Pr t « 0.7 is assumed, the 
calculations are in even better agreement with data. The reason 
could be due to a further reduction of turbulent mixing as a result 
of the wall being cooled. However, this effect has not been 
appropriately accounted for in the models, particularly their near- 
wall behavior. In other words, if highly cooled-waJl flows are 
to be predicted correctly, heat fluxes should be modeled 
separately rather than linking to momentum fluxes via a constant 
turbulent Prandtl number. 

fL 2 Discussion 

In the past, velocity profiles in wall coordinate were 
invariably plotted in terms of to illustrate the existence of the 
log-law and the constancy of k in compressible boundary-layer 
flows. The proposal was first suggested by van Driest [37] and 
later confirmed by Maise and McDonald [38] when they 
analysed ten sets of data in the Mach number range of 0 - 5. 

Since then, the compressible law of the wall is taken to be given 
by rather than by u + and k is considered to be about 0.41 and 
constant over the Mach number range of 0 - 5. The calculated 
and measured velocity plots given in Fig. 1 show support for the 
compressible law of the wall when it is written in terms of u + 
rather than u* Furthermore, k is determined to be 
approximately B.41 and is relatively constant over the Mach 
number range of 0 - 10. These results seem to conflict with the 
proposal of van Driest [37]. In order to resolve this seeming 
contradiction, the velocity plots of u7 versus In y* for cases 
55010504 (M« = 2.244, 0*/©, = 1) and 53011302 (M« = 
4.544, 0 w /0 r = 1) are shown in Fig. 13. In addition, the 
compressible law of the wall as given in Ref. 38 is shown for 
comparison. It can be seen that a line that is parallel to the 
compressible law of the wall can be drawn through a few of the 
data points spanning over a narrow range of y* . On the other 
hand, the calculated profiles are in agreement with data over a 
wider range of y* . The slopes of the calculated profiles are 
roughly parallel and are slightly larger than the slope of the 
compressible law of the wall shown. Therefore, irrespective of 
how the velocity profiles are plotted, the calculations are in good 
agreement with data. However, the slope of the log- law is given 
by 1/0.41 only when the profiles are plotted in terms of u + . 

The k and e equations for compressible flows are derived 
by assuming that there is no dynamical similarity between the 
compressible and incompressible fields. Therefore, the 
influences of fluctuating density on the mean and turbulence 
fields have to be accounted for in the modeled equations. This 
can be accomplished by first re-casting the exact k and e 
equations into forms that are similar to their incompressible 
counterparts. In other words, the viscous diffusion and 
dissipation functions have to be defined exactly like their 
incompressible terms. This procedure gives rise to additional 
terms in the k and e equations. These terms depend explicitly on 
compressibility and vanish when the fluid density becomes 
constant. One extra term in the k-equation is related to fluid 
dilatation and can be interpreted as compressible dissipation. 
The others are production terms that depend on the gradients of 
the mean pressure and mean viscous shears. All additional 
terms arc found to be relatively unimportant in the near- wall 
region, or 0 £ y* £ 50. This realization, therefore, allows the 
near-wall incompressible models to be extended directly to 
compressible flows without modifications, while still 
maintaining the balance of the modeled equations as a wall is 
approached. Models arc proposed for the additional terms in the 
k and e equations. The constants introduced by the new models 
are determined by calibrating the calculations against 
measurements in compressible flows. 

The near-wall two-equation model is used to calculate 
compressible flat plate boundary-layer flows with different wall 
thermal boundary conditions and free-stream Mach numbers. 



Comparisons are made with various mean flow measurements 
and with calculations of the k-co model. Good agreement is 
obtained between the present calculations and measurements. In 
particular, the log-law for compressible flows is recovered and 
the slope of the log-law is found to be fairly independent of free- 
stream Mach number for the range, 0 ^ M** £ 10, tested. Even 
though k-<n model gives a correct prediction of u + versus In y w , 
their velocity comparison in terms of <U>/Uoo versus y/5 shows 
substantial discrepancy with data. The discrepancy increases 
with increasing Mach number and can be attributed to a near- 
wall behavior that is not asymptotically correct. 

The following conclusions can also be drawn from the 
above analysis. Firstly, Morkovin s hypothesis is valid up to a 
free-stream Mach number of about 5 for flat plate boundary- 
layer flows with adiabatic wall boundary condition. This means 
that the effects of fluctuating density are becoming more and 

more important as M«, increases beyond 5. Secondly, the 
assumption of a constant turbulent Prandtl number is not 
appropriate for cooled wall thermal boundary condition. The 
reason is further reduction in turbulent mixing due to a cooled 
wall and this effect is not correctly accounted for in a constant 
turbulent Prandtl number approach. Most likely a heat flux 
model is required if the characteristics of cooled-wall 
compressible boundary-layer flows are to be predicted correctly. 
Thirdly, it is important to model the near-wall flow correctly if 
the overall boundary-layer characteristics are to be predicted with 
confidence. This point is substantiated by the k-e model 
calculations where all additional compressible terms in the 
turbulence equations are neglected. These results are in good 
agreement with measurements even though they differ slightly 
from the predictions of k-£ model/ 1 where all the additional 
terms are retained. In other words, an asymptotically consistent 
near- wall model is more important to the prediction of 
compressible boundary-layer flows than the inclusion of 
fluctuating density effects in the modeled equations. Fourthly, 
the predicted near-wall characteristics are very similar to those 
calculated for incompressible flows. In the range of free-stream 
Mach number tested, the calculated near-wall characteristics are 
essentially independent of Mach number and wall thermal 
boundary condition. Very near the wall, viscous diffusion of k 
is balanced by the dissipation of k. Beyond y* = 15, dissipation 
is balanced by mean shear production of k. In between these 
two regions, viscous and turbulent diffusion of k, production of 
k and dissipation of k are of importance in the budget of k. The 
additional compressible terms in the k-equation are essentially 
negligible in the near-wall region up to y* = 50. This is the 
reason why the model also performs well when the additional 
compressible terms are neglected in the equations. Finally, the 
term 3k/dy in the mean energy equation makes a significant 
contribution to the calculated Cf in the highly cooled wall case. 
Traditionally, this term is neglected. However, present analyses 
show that even though it is relatively unimportant in flows with 
adiabatic wall boundary condition, it cannot be neglected in 
flows with a highly cooled wall. The inclusion of this term 
degrades the prediction of C f . It is believed that the degradation 
is a result of an incorrect modeling of turbulent heat flux. 
Therefore, improvements should be directed at the relaxation of 
the constant turbulent Prandtl number assumption. 
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Case 

M. 

e^e r 

Model 


ayyXlO 4 

avexlO 7 

k + /e + (y*) 2 

55010504 

2.244 

1.0 

k-e model/ 1 

0.0987 

7.167 

-0.465 

0.50 

55010504 

2.244 

1.0 

k-e model/2 

0.0992 

7.198 

-0.998 

0.50 

53011302 

4.544 

1.0 

k-e model/1 

0.0824 

6.700 

-6.44 

0.50 

53011302 

4.544 

1.0 

k-e model/2 

0.0836 

6.760 

-11.79 

0.50 

73050504 

10.31 

1.0 

k-e model/1 

0.0741 

6.630 

89.5 

0.50 

73050504 

10.31 

1.0 

k-e model/2 

0.0771 

6.740 

-131.0 

0.51 

59020105 

5.29 

0.92 

k-e model/1 

0.0784 

6.120 

n.i 

0.50 

59020105 

5.29 

0.92 

k-e model/2 

0.0788 

6.140 

-5.88 

0.50 


Table 1. Asymptotic near- wall behavior of the turbulence properties. 











































































Figure 2. Comparison of calculated <U>/U*o with 
measurements. 


1 . Comparison of calculated u + with measurements. 









Figure 7. Near-wall distributions of Reynolds heat flux. 


Figure 8. Turbulent kinetic energy distributions across the 
boundary layer. 




Figure 11. Variation of the slope of k+ at the wall with 
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Figure 13 Mean velocity plots in terms of van Driest variable 
Figure 12. Skin friction coefficient variations with and + for cases 55010504 and 53011302. 










